MolDStruct: Modeling the dynamics and structure of matter exposed to ultrafast x-ray lasers with hybrid collisional-radiative/molecular dynamics.

We describe a method to compute photon-matter interaction and atomic dynamics with x-ray lasers using a hybrid code based on classical molecular dynamics and collisional-radiative calculations. The forces between the atoms are dynamically determined based on changes to their electronic occupations and the formation of a free electron cloud created from the irradiation of photons in the x-ray spectrum. The rapid transition from neutral solid matter to dense plasma phase allows the use of screened potentials, reducing the number of non-bonded interactions. In combination with parallelization through domain decomposition, the hybrid code handles large-scale molecular dynamics and ionization. This method is applicable for large enough samples (solids, liquids, proteins, viruses, atomic clusters, and crystals) that, when exposed to an x-ray laser pulse, turn into a plasma in the first few femtoseconds of the interaction. We present four examples demonstrating the applicability of the method. We investigate the non-thermal heating and scattering of bulk water and damage-induced dynamics of a protein crystal using an x-ray pump-probe scheme. In both cases, we compare to the experimental data. For single particle imaging, we simulate the ultrafast dynamics of a methane cluster exposed to a femtosecond x-ray laser. In the context of coherent diffractive imaging, we study the fragmentation as given by an x-ray pump-probe setup to understand the evolution of radiation damage in the time range of hundreds of femtoseconds.

Shell-dependent photofragmentation dynamics of a heavy-atom-containing bifunctional nitroimidazole radiosensitizer.

Radiation therapy uses ionizing radiation to break chemical bonds in cancer cells, thereby causing DNA damage and leading to cell death. The therapeutic effectiveness can be further increased by making the tumor cells more sensitive to radiation. Here, we investigate the role of the initial halogen atom core hole on the photofragmentation dynamics of 2-bromo-5-iodo-4-nitroimidazole, a potential bifunctional radiosensitizer. Bromine and iodine atoms were included in the molecule to increase the photoionization cross-section of the radiosensitizer at higher photon energies. The fragmentation dynamics of the molecule was studied experimentally in the gas phase using photoelectron-photoion-photoion coincidence spectroscopy and computationally using Born-Oppenheimer molecular dynamics. We observed significant changes between shallow core (I 4d, Br 3d) and deep core (I 3d) ionization in fragment formation and their kinetic energies. Despite the fact, that the ions ejected after deep core ionization have higher kinetic energies, we show that in a cellular environment, the ion spread is not much larger, keeping the damage well-localized.